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Abstract 

We present transport measurements of tunnel junctions made between Cu and graphene in a 
magnetic field. We observe a transition to a Landau level like structure at high fields, as well 
as a set of sharp features in the tunneling spectra that shift with gate and tunnel probe voltage 
along the lines of constant charge density. We explain the sharp features with the formation of 
degeneracy split localized Landau levels, and addition of electrons to those levels one by one. A 
large capacitive coupling to the tunnel probe also increases the gate voltage spacing between the 
Landau levels. 



The electron and hole bands of pristine graphene exhibit linear dispersion, and meet at 
a point in reciprocal space of zero energy and vanishing density of states (DOS) known as 
the Dirac point. This situation leads to charge carriers that behave as massless relativistic 
particles, and a large number of theoretical and experimental studies have investigated the 
consequences of this fact. For review see Refs. 1,2 . 

Of increasing importance for future graphene technology is the interaction between 
graphene and a metal. Since metals will eventually be involved in some way in graphene 
circuits, as interconnects or gates, it is critical to understand their effect on the graphene. 
Theoretical 3 ' 4 and experimental 5 " 10 publications have found that one of the major effects of 
a metal on graphene is doping. In addition to this, it is believed that some metals interact 
weakly, leaving the linear dispersion of graphene undisturbed, and thus it may be desirable 
to use such metals near regions where graphene properties are most important. 

In this work, Cu/graphene tunnel junctions are used to investigate the structure of the 
field dependent graphene density of states (DOS) beneath a metal contact near the Dirac 
point. We can also identify certain aspects of the data which can only be explained by 
the more complex physics of single electron charging, providing physical evidence for the 
strongly localized Landau levels which play a crucial role in the quantum Hall effect. 

Previous studies of tunneling in graphene in magnetic field have been performed with 
Scanning Tunneling Spectroscopy (STS) 11-13 , however these studies cannot study graphene 
in proximity to a bulk metal. STS has also been performed on graphene on metal substrates 
such as Cu 14 and Au 15 , but the lack of a gate electrode or magnetic field in these studies 
creates difficulty in interpreting peaks in the dl/dV spectrum. Finally, solid state devices in 
a geometry similar to our own have been performed with a superconducting Pb electrode. 16 
This study also claims, as we will later, that a localized region is created underneath the 
metal electrode. However, the study investigates excitations within the superconducting gap 
region in zero field, whereas our present study concerns tunneling between a normal metal 
and graphene in the presence of up to 12T. 

In order to fabricate our tunneling junctions, graphene is first mechanically exfoliated 
from single crystal graphite obtained from Naturally Graphite™ onto a highly doped Si 
wafer with a 300nm layer of thermally grown oxide, 17 and then single layers are confirmed 
by Raman spectroscopy 18 . Cu junctions were directly evaporated onto graphene by e-beam 
lithography. After lift-off in acetone, the device is placed in dilute nitric acid (~12% HN0 3 ) 
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for several seconds, and then rinsed with DI water, isopropanol, and blown dry. The acid 
promotes oxidation of the Cu, and some junctions become highly resistive after this treat- 
ment (>1 MQ), afterwards they continue to age, though slowly, taking several weeks to 
fully oxidize. At low temperatures, the resistance values are stable, as the oxidation rate is 
greatly suppressed. See Ref. 9 for additional details about fabrication. 

Narrow junctions (< 200 nm) were used to insure that the acid could partially penetrate 
the Cu/graphene interface before the lead was overetched. Though slight tarnishing could 
be seen in the leads, only leads that were highly conductive at low temperatures (kfis) 
were used as a grounding electrode. The conductance of the tunnel probe was not severely 
suppressed at low temperatures, meaning that the junction between metal and graphene 
provided the primary mode of resistance in the circuit, and not a discontinuity in the metal 
lead. Also, as we shall show, the tunneling spectra exhibit expected graphene properties. 

We report on one junction in particular, however the primary results presented here were 
observed in two other devices that displayed tunneling properties. Fig. 1-A shows an optical 
image of the completed device, with a dotted line encircling the particular probe used in 
data discussed. Many of the other probes on this device did not become highly resistive after 
the short exposure to HN0 3 . In the present case the junction was only partially overlapping 
the sample, and in two other cases the junction was entirely overlapping. All measurements 
presented here were taken while the sample was in vacuum in a cryostat operated at 2.2 K, 
and equipped with a superconducting magnet. Tunnel spectra were gathered by applying 
a DC probe voltage, V p , across a highly resistive and an ohmic probe, and a gate voltage, 
V g , between the Si back gate and an ohmic probe, and measuring the current as shown in 
Fig. 1-B. A numerical derivative was then taken to calculate the conductance. 

Our data consists of mapping the conductance versus V g and V p at different magnetic 
fields. Fig. 2 displays our conductance maps at 0, 6, 9, and 12T. The tunnel spectra 
are dominated by the gate voltage independent zero bias anomaly (ZBA), similar to other 
tunneling experiments. 12 ' 19 We can analyze the data despite the ZBA by inspecting how 
tunnel resonances in the spectra shift with V g . Those features that shift with gate voltage 
are related to the graphene DOS at energy eV p + EF(V p ,V g ), and their slope is directly 
proportional to the DOS, 9 where E F (V p: V g ) is the Fermi level in graphene. By applying a 
gate or probe voltage, charge is added to or subtracted from the graphene. In regions of 
high DOS, the added charge does not shift the Fermi level much, whereas in regions of low 
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DOS, the added charge shifts the Fermi level quickly. 

Looking at Fig. 2-A, lines of high conductance can be seen at OT. These are resonances in 
the DOS caused by the disorder potential, and the slope dV g /dV p becomes slightly flatter at 
positive gate voltage, where the density of states is lower. This we interpret as a broadened 
Dirac point, 9 and its location at positive gate voltage means the graphene is hole doped. 
As the magnetic field is increased, these lines are seen to break up, and are replaced by 
lines with a staircase like structure of alternating high and low sloped lines. The Staircase 
like structure is indicated by the guides in Fig. 2-E and F. We identify these lines as the 
formation of Landau levels. Since our junction covers a large area compared to the expected 
magnetic length in graphene, as well as the charge puddles caused by inhomogeneous doping, 
it is possible that the conductance maps represent conductance from multiple channels, and 
so features in addition to those predicted by theory may be present. 

To compare our conductance maps with those expected from perfect graphene we nu- 
merically calculate the tunneling conductance assuming constant transmission with respect 
to energy for states between the Fermi level and the probe voltage for given magnetic field, 
probe, and gate voltages. The semiclassical model for the ideal Landau level structure of 
graphene gives the DOS in a magnetic field as X^-oo ^tt^ {p ~~ vpsign(i)^j2ehB\i\^ . The 
use of such a DOS is justified by analytic 20 and numeric 21 calculations of the electronic spec- 
trum for confined graphene quantum dots in a magnetic field, which show that the localized 
states of the quantum dot still converge upon the expected landau level spectrum at high 
magnetic fields. 

To convert from the DOS to the conductance map, we use the capacitor model described 
in Ref. 9 to solve for EpiVg, V p ). We arrive at Fig. 3-C by applying a 3.5meV Lorentzian 
broadening function to the ideal DOS and using a capacitance ratio between the tunnel 
probe and the back gate's capacitance per unit area of 250. A staircase structure can be 
seen to emerge in which the gate voltage is unable to shift features in the spectra until the 
highly degenerate Landau levels are filled, and then the spectra are quickly shifted when the 
DOS is low. 

We can use the same assumptions to simulate conductance maps vs. B and V p . Ideally, 
such a map should show a series of parabolas shifted from the center depending on the 
doping level and applied gate voltage. The simulation of a conductance map vs. B and V p 
in Fig. 3-A shows that the expected parabolic lines are disrupted by shifts in the Fermi level 
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caused by the tunnel probe. We show for comparison data taken from our tunnel probe 
which shows similar behavior (Fig. 3B), at V g = OV. 

An important observation of this paper is the sharp resonances that begin to appear in 
the upper portion of the 6T map as seen in Fig. 2-B, and become very well resolved, with the 
lowest resonance appearing at V g = —20V and V p = OmV at 12T in Fig. 2-D. In prior work, 9 
we showed that the low carrier density regions in graphene, produce anomalies in the tunnel 
conductance maps that shift along the lines of constant charge density C p V p + C g V g = const. 
The resonance shifts in Fig. 2B-D can be attributed to low density regions with C p /C g = 250. 

We interpret these resonances to be due to single electron charging effects on localized 
states near the filling of Landau levels. In the progression from Fig. 2-B-D, the set of sharp 
resonances where the lowest resonance appears near +25 V do not appear to shift, whereas 
the second set of resonances is seen to shift from near 0V in Fig. 2-C to ~-20V in Fig. 2-D. 
This, combined with its clear appearance before any other level, leads us to the conclusion 
that the Dirac point is near +25V, and that these resonances are associated with the n=0 
Landau level. The other set of resonances which show near -20V at 12T would therefore 
be the n= -1 Landau level. The bright, high conductance, positive sloping lines at -20V, 
indicated by white guides to the eye in Fig. 2-F, changes slope after intersecting with the 
top most resonance line. The change in slope is consistent with that found in Ref. 12 , which 
is another indicator that the resonances form when the LL level is filling. 

The gate voltage spacing is much larger than expected if these are in fact the n=0 and 
n=-l levels. The spacing between the bottom resonance of the n=-l, and the bottom 
resonance of the n=0 is ~45 V, when 16 V is expected at 12T. A possible explanation is 
a capacitive division due to the tunnel probe, and numerical solution of the electrostatic 
equations shown in Fig. 3-C demonstrate that the spacing in gate voltage between the 
various Landau Levels is nonuniform. Since graphene has a low density of states compared 
to that in a metal, the electric potential of the graphene can shift without a corresponding 
change in the electrochemical potential of the grounding electrode as long as the Fermi level 
of the graphene compensates such that fiF + e<p = 0. In ideal graphene in a magnetic field, 
the density of states between the Landau levels is zero, and the electric potential follows 
capacitive division between the gate and the probe potential. 

The required gate voltage to cause a given shift in graphene's electric potential can be 
solved using the equation C g (AV g — A0) + C p (— A0) = — AQ where AQ is the induced 
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charge on graphene, A0 is the electric potential shift in graphene under the probe, and the 
tunnel probe potential is held at zero. To shift from the n=-l to n=0 Landau level requires 
first that the electric potential is shifted such that the Landau level to be filled aligns with 
the electrochemical potential of the leads meaning that — eA0 = A/ip = vpVZhB = AeiL, 
and second that the necessary charge density to fill a Landau level, —4e 2 B/h, is supplied 
from the leads. Therefore, eAV g = ^s+^a) Ae LL + ^j^-- The first term on the right is small 
when C p is small (~ l25mV at 12T) and can typically be ignored, but if C p /C g = 250 then 
this term becomes 31V, and must be taken into account. Thus AV g = 47V as in our data, 
rather than the expected 16V. The resonances near 0V in Fig. 2-C are not well developed 
enough to accurately measure their position. 

A closer inspection of Fig. 3-D reveals the resonances are spaced by approximately lOmeV 
along the V p direction. The fact that the spacing is much larger than the expected Zeeman 
(1.4meV) splitting in the 1st Landau level of graphene, indicates that the charging energy of 
the Landau level dominates these spacings. It vaguely appears that the tunnel conductance 
of the sample, near zero probe voltage, is suppressed within the region enclosed by four dia- 
monds, indicated by the arrows in Fig. 3-D, suggesting that the resonances can be explained 
by the theory of Coulomb blockade. 22 In the present case, a localized state within a Landau 
level becomes the central island, tunnel coupled to both the metal lead and an extended 
state within the graphene which in turn couples to the grounding electrode. 12 Such samples 
exhibit sudden onsets of current along the sides of a Coulomb diamond. The two negative 
slopes of the Coulomb diamond (AV p /AV g = —C p /C g ) occur when the island is at a large 
enough potential for an electron to tunnel into (out of) the extended state, followed by an 
electron tunneling off (onto) the probe. The spacing between the negative sloping lines, 
measured parallel to the probe voltage axis, is e/AC p , where A is the area of the island, 
assuming the energy level spacing due to Zeeman and the valley degeneracy splitting much 
smaller than the charging energy. 

We can derive the same expression for the spacing more generally, without using the 
theory of Coulomb blockade, which is important because the Coulomb diamonds are rather 
vague. As discussed earlier, the regions of low density in the conductance maps versus C p 
and Cg, shift along the lines of constant charge density q = C P V P + C g V g . According to our 
interpretation, moving from one line to the next one, along the direction of increasing V p 
(V g ), corresponds to adding an electron to an unoccupied, degeneracy split LL. Thus, q to 
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changes by e/A from line to line, leading to the spacing e/AC p (e/AC g ). 

Thus, we obtain AC V = 16aF. Assuming a gate capacitance per unit area of 124/iF/m 2 
as measured on a test wafer from the same batch of oxidized Si wafers, as well as C p /C g as 
described earlier, we calculate a radius of the localized state of 13 nm. The magnetic length 
of graphene at 12T is 7.4nm 23 , which suggests that tunneling could indeed be occurring 
through a single localized Landau level. We note that the resonances located near +25V 
gate voltage have a larger spacing (~15meV). As the region of graphene affected by the gate 
becomes much smaller than the 300nm gate dielectric thickness, the breakdown of the plane 
capacitor model can explain small discrepancies in the measured and expected length scales. 

Our data is to be compared with, 12 in that work, the Landau levels are identified as a 
staircase structure in the conductance maps, and at high fields, sharp resonances are seen to 
emerge at the filling and emptying of Landau Levels. It was suggested that the larger energy 
spacing near the zero Landau level was due to a lifted valley degeneracy, adding lOmeV at 
8T to the charging energy in the n=0 Landau level. The staircase structure has recently 
also been observed in ref. 13 . 

In conclusion, transport through graphene tunnel junctions in a high magnetic field 
demonstrates localized states in the quantum Hall regime. The direct probing of a localized 
state in a Landau level is both important as confirmation of earlier experiments, and a 
new way to study Landau level physics in simple solid state device. Also, the capacitive 
coupling between the Landau level and the probe electrode increases the apparent spacing 
between the Landau levels in gate voltage, and should be accounted for when assigning 
energy spacings to Landau levels in spectroscopic studies. 
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FIG. 1. A: An optical picture of the graphene used in this experiment, the tunnel probe is circled 
with a dotted line. B: A cartoon of the experimental setup showing the tunnel probe with a 
localized and extended state beneath it. 



10 



10x10" 9 20 30 40 




-20x1 0" 3 20 -20x1 0" 3 20 

Probe Voltage Probe Voltage 



FIG. 2. A-D: Progression of conductance maps from to 12T. The original set of resonances, due 
to the disorder potential, are broken up and replaced by a series of staircases. E-F: The same 
conductance maps as C & D, with white guidelines to mark the formation of staircase structure 
in the resonances, and black lines to show the formation of structure associated with the charging 
of localized states. The white line in F decreases in slope exactly at the point where it crosses the 
black line showing that the DOS drops after filling the Landau level, consistent with 12 
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FIG. 3. A: A simulation of a V p vs B map assuming undoped graphene and -35 V on the back gate. 
B: Data taken at 0V gate voltage in the sample under discussion. C: A simulation of a V p vs V g 
map, showing the staircase structure of the Landau levels at 12T, a box denotes the approximate 
range of at least one of the conductance channels measured in our data. D: Arrows mark the 
Coulomb diamonds at 12T. 
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